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Primary lymphoid follicles in secondary lymphoid tissue of mammals are the backbone 
for the formation of follicular dendritic cell networks. These are important for germinal 
center reactions. In the context of organogenesis molecular requirements for the formation 
of follicles have been identified. The present study complements this work with a simulation 
of the dynamics of the primary lymphoid follicle formation. In contrast to other problems 
of pattern formation, here, the homeostasis of the cell population is not governed by a 
growth-death balance but by a flow equilibrium of migrating cells. The influx of cells 
into secondary lymphoid tissue was extensively studied while less information is available 
about the efflux of lymphocytes from secondary lymphoid tissues. This study formulates 
the minimal requirements for cell efflux that guarantee a flow equilibrium and, thus, a 
stable primary lymphoid follicle. The model predicts that in addition to already identified 
mechanisms a negative regulation of the generation of follicular dendritic cells is required. 
Furthermore, a comparison with data concerning the microanatomy of secondary lymphoid 
tissues yields the conclusion that dynamical changes during the formation of FDC networks 
of the lymphatic endothelium are necessary to understand the genesis and maintenance of 
follicles. 
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1 Introduction 



The primary lymphoid follicle (PLF) is the precursor state of germinal centers in secondary lymphoid 
tissue (SLT) [1]. The most important constituents are naive B cells and follicular dendritic cells 
(FDC). PLF serve as a filter for antigen and bring B cells into contact with antigen presented by FDC 
in order to start a germinal center reaction [2]. 

In mammals PLF first develop around birth [3—8] . During the past decade many molecules relevant 
to the formation of PLF could be isolated. With an extensive number of knockout and antibody- 
blocking experiments the contributions of TNF-a, LTai/?2, LT«3, LIGHT, and TRANCE secondary 
lymphoid organ structure and in the formation of PLF have been identified (reviewed in [9-11]). 

The homeostatic chemokines CCL19, CCL21 and CXCL13 are involved in guiding the lymphocytes 
to their compartment (reviewed in [12]). CXCL13 is secreted by FDC and acts as chemoattractant 
for B cells via the CXCR5 chemokine receptor. CCL19 and CCL21 are produced by different cell 
populations in the T zone of SLT. The common receptor CCR7 is expressed by T cells. The expression 
of these chemokines and the presence of FDC depend on LTai/?2, TNF-a and related molecules. 
The precise molecular interactions and chemokine expression profiles slightly vary for different SLT 
[9—11], however, in most SLT a persistent LTai /^-stimulus is needed for both the induction and 
the maintenance of FDC networks [13]. Further studies suggested that the maintenance of the PLF 
structure is mediated by a positive feedback loop [14]: B cells are stimulated by CXCL13 to express 
high levels of LTai/32- This in turn stimulates FDC to produce CXCL13. 

An interesting property of PLF is their equal size of several hundred /im in all mammals ranging 
from mice to horses [15-22]. Thus, the follicle size does not scale with the size of the animals or of 
the organs under consideration, which is considered to be a strong boundary condition for simulations 
of PLF formation and maintenance. The formation of PLF is studied with a simulation based on 
a previously introduced agent-based model [23] on top of a regular triangulation [24-26]. A clear 
separation of B and T cells is achieved using chemoattraction to an exit spot for the lymphocytes 
overlaying the different chemokine responses of B and T cells to their homing chemokines CXCL13, 
CCL21, and CCL19. Formation of a PLF and an adjacent T zone is found when assuming the 
generation of FDC from T zone stromal cells by sufficiently large B cell aggregates. However, then, 
the location of the PLF relative to lymphatic vessels is not in agreement with experiment. The present 
study concentrates on the correct formation of the lymphatic vessels and the PLF without considering 
the dynamics of the T zone of T cells. The model includes lymphangiogenesis, which dramatically 
changes the dynamics of the PLF formation, its geometry and also the regulation mechanisms involved 
in PLF maintenance and size regulation. It is found that a negative regulation of the FDC generating 
process is required. Furthermore it is demonstrated that a chemotactic activity of B cells for SIP is in 
contradiction to PLF morphology, suggesting that SIP chemotaxis is not active during PLF formation 
and maintenance in vivo, which was also found in experiment [27]. 

2 Model of primary lymphoid follicle formation 

This section describes the properties of a minimal model for PLF formation. In addition the model 
of lymphangiogenesis is discussed. The mathematical framework is briefly summarized in the method 
section (Sec. [5]) and is presented in more detail in [23]. 
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2.1 Flux of naive B cells 



Essential ingredients to study the formation of the PLF are the B cell flux, the FDC and the FDC 
precursor cells. The B cells are constantly entering SLT via high endothelial venules (HEV) in the cases 
of mucosa-associated lymphoid tissue (MALT) and lymph nodes (LN) or along the central arterioles 
through the peri-arteriolar sheath of the spleen (reviewed in [12]). The exit route of lymphocytes 
through lymphatic endothelium (LE) is less clearly identified. Recent experimental data (see [12]) 
suggest that B cells leave MALT via efferent lymphatic vessels (ELV) . In LN and spleen the lymphatic 
sinuses guide lymphocytes to draining ELV. As the sinuses merge with the ELV they are considered 
to belong to the efferent lymphatic system in the simulation. 

The morphology of the LE is very similar across the different SLT and species. Studies using 
corrosion casts revealed that the LE is formed around the follicles. The precise shape ranges from 
half-open baskets around the follicle base to almost closed shells enclosing the whole follicle. The LE 
forms a dense network in the non-follicular area [21,28-37]. 

The role of LE in lymphocyte emigration is supported by a study of the corresponding structure 
in the Bursa of chicken [38] and the identification of openings in the LE nearby follicles [28,31]. The 
exit of lymphocytes via LE is further supported by a newly identified mannose receptor on LE. This 
receptor may guide the specific exit route of lymphocytes via the lymphatic sinuses [39]. It is almost 
exclusively found on the lymphatic vessels in the LN and binds CD62L, thus, guiding lymphocytes to 
the ELV and allowing their exit from LN. 

SIP might be involved in the regulation of lymphocyte egress from SLT [12,40-42]. SIP acts as 
a chemoattractant for lymphocytes in vitro [40]. Thus, the dynamics of the SIP receptor SIPi is 
considered in the simulation: The LE is widely engulfing follicles such that lymphocytes have to 
migrate around or across these structures in order to reach the follicles. Therefore a minimal transit 
time through SLT should be governed by some mechanism that prevents the lymphocytes from entering 
and transmigrating across the LE. In the simulation it is assumed that SIPi is downregulated when 
B cell enter SLT and required for lymphocyte emigration. This is supported by studies using the 
SIPi agonist FTY720 [12,40-42]. There is good agreement between the minimal transit time of 
lymphocytes [33,43-45] and the time required to fully upregulate SIPi on lymphocytes [42] both of 
which are in the order of 3 hours. 

2.2 Origin of FDC 

One fundamental question in the PLF formation is: Where do the FDC come from? The most 
commonly accepted view is that FDC are derived from stromal cells which might be related to the 
stromal cells observed in the T zone. These stromal cells are fibroblastic reticular cells (FRC) [46-49]. 
Comparative reviews supporting the mesenchymal origin of FDC can be found in [50]. The evidence 
for the relation of FDC to FRC is given by shared markers [51]. There exists a gradual decrease 
of the fibroblast marker ASO-2 from the T zone to the heart of a germinal center. This decrease is 
accompanied by a gradual morphological change from the 'classical fibroblastic' to the typical FDC 
morphology [52]. A similar change of marker expression has been observed in the PLF although a few 
intermediate steps are missing. Indirect evidence for the reverse transition from FDC to stromal cells 
is provided by culture experiments in which purified FDC gradually loose characteristic FDC marker 
one by one, finally turning stromal cell-like again [53-55]. This suggests that persistent stimuli are 
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required to maintain the FDC phenotype. Thus, the study of typical FDC marker expression indicates 
a smooth transition from stromal cell precursors to FDC by subsequent accumulation of 'FDC-ness' 
of the stromal cells. A recent study of FDC development in human LN supports the idea of local FDC 
differentiation: FDC seem to acquire more and more markers like those found in culture experiments 
extending the 'FDC-ness' picture to several other markers [20]). 

Observations of slow replacement of FDC by migratory precursors and of inefficient seeding by bone 
marrow-derived cells argue against a hematopoietic origin of FDC [56]. However, a clear relationship 
between FRC and FDC could not be identified in a study based on newly developed markers [49]. 
Bone marrow transfer into severe combined immunodeficiency (SCID) mice revealed donor-derived 
FDC suggesting that the stromal cells may originate from the bone marrow [57]. TNFR-I knockouts 
and mice given injections of a blocking TNFR-FIg exhibit immature FDC or FDC precursors around 
splenic follicle suggesting migrative FDC progenitors [58]. Similar conclusions have been drawn from 
TNF-a~/~ and TNFR-I - / - [59]. More support for the immigration of FDC precursor instead of local 
differentiation comes from the expression pattern of TuJl, a microtubule-protein [60]. The broad 
distribution with single positive cells scattered around lymphoid follicles suggests migratory activity 
of FDC. Reconciling both theories, FDC might be a heterogeneous population that can develop from 
hematopoietic and stromal precursors [56,61]. This favor the vantage point adopted in this study 
considering FDC as a state of cells rather than a differentiated form of some cell type. 

In the simulation FDC differentiate from FRC. This can be viewed as FDC being an 'excited' state 
of stromal cells. However, this must not be interpreted in a strict sense. The mechanism in the 
simulation is indistinguishable of the interpretation that FRC and FDC have a common progenitor 
which can develop in either FRC or FDC depending on the external stimuli. Newly developing FDC 
replace existing FRC at their site of generation. 

In a simplified PLF model it could be shown that sufficiently large aggregates of B cells forming 
around stromal cells can induce the presence of FDC if basal levels of LTai/?2 are expressed on the B 
cells [23]. The migration properties of B cells derived from the experiment are such that the duration 
and size of the B cell-FDC contacts are sufficient to allow for a high LTai/?2 stimulus of a single 
FDC by the interaction with many B cells. The stimulus can last for several hours which seems to 
be reasonable according to the experimental data [13,23,29,62]. Thus, we conclude, that the positive 
feedback loop [14] (which is included in the simulation) is not absolutely required to explain initial 
FDC formation. 

2.3 Phenomenological lymphangiogenesis model 

Lymphangiogenesis is included in the model in order to generate a LE distribution that matches the 
experimental situation [21,28-37]. It is assumed that the LE is not preformed which would determine 
shape and position of PLF a priori. A preformed LE gap would immediately raise the question what 
generates this specific LE pattern. Thus, in the simulation the generation of FDC is anti-correlated 
with the LE dynamics: LE is degraded when FDC are formed in the vicinity. It is assumed that the 
anti-correlation between FDC and LE is local, i.e. the LE-degradation process is not governed by a 
long ranged diffusion of molecular messengers. In a similar manner LE is forming at places where 
FDC have transformed back into FRC. Angiogenesis and lymphangiogenesis models are still under 
development [63-66]. Therefore, the present tissue simulation takes into account vessel formation and 
degradation in the described phenomenological sense without considering the underlying mechanisms. 
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2.4 Internalization of CXCR5 



The B cells have to leave the PLF before they can exit from the SLT. LE will dynamically enclose the 
follicles. Thus, B cells frequently reach the surface of the follicle and randomly move until a lymphatic 
vessel has been found for exit. Alternatively, they might return to the follicle. It is assumed that this 
search for exit points is regulated by B cell chemotaxis. As the chemotactic response of leukocytes 
can be modulated by receptor internalization (reviewed in [67]), in the simulations this regulation is 
mediated by internalization of the CXCL13 receptor CXCR5. 

2.4.1 Regulation of chemotaxis 

Studies on neutrophils demonstrated that desensitization of chemoattractant receptors occurs predom- 
inantly for high concentrations of the chemoattractant [68] . One mode of desensitization is the inter- 
nalization of the chemokine receptor. However, in line with the dependence on high concentrations, it 
has been shown that internalization is not absolutely required for chemotactic responses [67,69-71]. 
For example CCL21 does not induce the internalization of CCR7 in contrast to CCL19 which causes 
the downregulation of this chemokine receptor [70]. It is a general observation that different ligands 
of the same receptor cause different internalization levels [67]. 

Of note, some experiments fail to detect CXCL13 responses of naive B cells freshly isolated from 
tonsils although the B cells where equipped with high levels of the corresponding receptor CXCR5 [72]. 
Thus, the presence of chemokine receptors CXCR5 is necessary but not sufficient to cause chemotaxis 
of naive B cells. Most likely the suppression of the function of the CXCR5 chemokine receptor is 
mediated by a desensitization mechanism either by internalization or cross-desensitization (see below). 

The importance of the receptor state can be demonstrated by the response of B cells towards 
CCL21. The surface CCR7 level on freshly isolated B cells is undetectable [73]. The protein can be 
identified in the cytoplasm and is brought to the surface in chemokine-free culture. This suggests 
receptor internalization by natural exposure to CCL21 or CCL19 in the tissue. Freshly isolated T 
cells also show little attraction by CCL19 [74] suggesting an internalization mechanism for T cells as 
well. However, the CCR7 levels have not been investigated in detail. Internalization may explain why 
despite the presence of CCL21 in the T zone, and lymphocyte expression of the associated receptor 
CCR7, no chemotaxis has been observed in this area [75]. 

The presented data supports the notion that the chemotactic response of lymphocytes in SLT can 
be modified by receptor internalization. Therefore cells may be unresponsive to chemokines despite 
proper receptor expression and in vitro responses. 

Note, that within the model framework regulation of chemotaxis on the basis of other mechanisms 
like transcriptional regulation is not incorporated. 

2.4.2 Complex responses to multiple chemokines 

In several systems cross-talks between different chemokines/chemokine receptors have been observed 
[71]. An example is the cross-desensitization that has been reported for CXCL12 which blocks the 
CCL19 response of lymphocytes while CCL19 does not block CXCL12 [76]. As expected cross- 
desensitization is found when different chemokines use the same receptor - as shown for CCL19 
and CCL21 [77]. 
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Not in all cases the effect of multiple chemoattractant is pairwise blocking of the responses. It is 
also possible that the cell computes a vector sum of the incoming signals to determine an average 
direction of multiple chemoattractants [78]. Only the direction of cell migration is influenced by the 
integrated response, while the speed remains unchanged. This reaction is altered if multiple signals 
are given in some sequential order. Then the direction is dominated by the newest chemoattractant 
even for lower concentrations and/or gradients [78]. An explanation for this behavior would be the 
desensitization of receptors for instance by internalization (see above and [70]). The desensitization 
state represents some memory for the chemoattractants a leukocyte is encountering. 

Instead of desensitization a hierarchy of chemoattractants might exist. Indeed, two signaling path- 
ways for two different chemoattractant receptors have been identified [79]. The hierarchy levels are 
named 'target' and 'intermediary' chemoattractants. When only concurrent intermediary signals are 
provided the cell responds by computing the vector sum otherwise the target chemoattractants are 
preferred. 

Which type of response results from multiple signals may rely on the signaling cascade, i.e. which 
receptors share or use a concurring cascade to induce directed cell migration. The data available on 
the chemokines CCL19, CCL21, CXCL13, and SIP, which are important for the PLF system, are 
not sufficiently conclusive to allow the determination of a signaling hierarchy. Activated naive B cells 
upregulate there CCR7-levels by a factor of two to three and consequently enhance their response to 
CCL19/CCL21. This balances their CXCL13 response such that they migrate to the border between 
the PLF and T zone [10,12,72,73,80-82]. When either CXCR5 is overexpressed or CCR7 is lacking 
B cells fail to relocate to the follicle border upon antigen-stimulation and remain in the follicle [81]. 
Similar if CCR7 is overexpressed by genetic manipulations the B cells locate at the border without 
antigen and are moving farther into the T zone. Analogously, activated T cells upregulate the response 
to CXCL13 relocating to the follicle border as well [12,80,83,84]. 

Overall this suggests that the relevant chemokines in the PLF system have the same hierarchy 
level and the chemotactic response of lymphocytes is a vector sum of concurrent chemokine signals. 
Therefore this mode of response for CCL19, CCL21, CXCL13 and SIP has been chosen in the model: 
The speed of lymphocytes is set to a constant value derived from two-photon imaging experiments 
[75,82,85,86]. The direction of the chemotactic response is the weighted average direction of the 
chemokine gradients. The weight is provided by the gradient of the bound chemokine at both ends of 
the cell assuming that the cell senses the difference of bound molecules across its diameter. 

3 Results 

Simulations are performed on the basis of the model as introduced in the previous section. The 
minimal set of assumptions entering the model are 

• fixed number of non-migrating stromal cells, 

• constant entrance of B cells via small number of HEV, 

• B cells use the LE for egress from the SLT, 

• generation of FDC by B cell LTai/?2 signaling, 
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Figure 1: Three dimensional slice projections from simulations (rendered with POVRay, 
http: / / www.povray.org[ ) . Objects are shown as colored spheres (yellow: FDC, green: FRC, 
white: naive B cells, red: HEV, dark grey: ELV). The internalization process destabilizes 
the follicle shape, even when the FDC network is of fixed size and ELV are distributed in the 
T zone (A-D, interval between images is 1 h). Also the dynamics of the efferent lymphatic 
destabilizes the follicle shape (E: 5 h, F: 122 h, G: 5 h, H: 57 h). Although the morphology 
of the follicle is correct with the LE surrounding a spherical follicle the situation does not 
equilibrate. This is independent of the density of the LE (E,F: low density; G,H: high den- 
sity). Note that at high LE density the distinct spheres representing the exit spots merge 
into a seemingly homogeneous background in the pictures. Introducing negative regulation 
with free FDC dynamics the follicle forms around the HEV (I: 10 h, J: 60 h). If however, like 
in the shown case, B cells chemotax towards CCL21 the follicle may elongate towards the 
center of the stromal network (J). Simulating B cell homeostasis with negative regulation 
of FDC generation and a preformed FDC network with high density LE the follicle rapidly 
adopts its final shape and keeps its position and does not relocate to the HEV (K: 5 days). 
When SIP is acting in the dynamic PLF formation a disturbed follicular border occurs (L: 
1 day). Also the PLF appears large due to the lower density of B cells while the number of 
B cells is not significantly increased (not shown). 
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replacement of FRC by generated FDC, 



• removal of FDC (replaced by FRC) in case of lack of LT«i/32 stimulation, 

• secretion of CXCL13 by FDC, 

• secretion of CCL21 by FRC, 

• secretion of SIP by LE, 

• chemotactic activity of B cells to CXCL13, CCL21, and SIP, 

• B cells leaving the SLT via LE when SIPi levels are sufficient to dominate CXCL13 response. 

These assumptions alone lead to reasonable cell dynamics. Restricting the LE to a small area generates 
a PLF that is stable in size, shape, and position (see [23]) using physiological values for the parameters 
(Table[T]). However, the position of the follicle is in contradiction to microanatomical data [21,28-37]. 
Introducing the correct relative position of PLF and LE with dynamic LE, dramatically changes the 
follicle stability. Thus, the stability of follicle shape and size in an equilibrium of cell flow crucially 
depends on further mechanisms. 

In the following possible mechanisms are discussed that guarantee the stability of the follicle. The 
most intuitive hypothesis, is related to internalization dynamics of chemotaxis receptors. It is shown 
that this hypothesis leads to unphysiological results and is therefore unlikely to be the relevant factor 
of physiological follicle formation and maintenance. Realistic results could only be achieved by the 
assumption of a so far unknown mechanism of negative regulation of FDC generation. This mechanism 
is a prediction of the model. 

3.1 Including internalization dynamics 

The internalization of CXCR5 (modeled by Eq. [5] and [23]) destabilizes the shape of the follicle defined 
by the B cells (Fig. [TJ A-D). In a first step, the dynamics of the FDC network was switched off in 
order to isolate the effect of B cell internalization dynamics. 

The instability is caused by uptake of CXCL13 when the CXCR5-CXCL13 complex is internalized: 
The chemotactically responding B cells act as a sink for CXCL13. This generates steep gradients in the 
chemokine distribution which guide the B cell movement. Interestingly, these local gradients, induced 
by the responding lymphocyte population itself, can become so strong that they reverse the chemotaxis 
gradient. Then cell movement is not strictly directed towards the chemokine source anymore. 

In principle, a large diffusion constant could counter-balance this tendency. However, the reversal of 
chemokine gradients is most pronounced at larger distances of the chemokine source. An unphysiolog- 
ically large diffusion constant would be required to counter-act the local dynamics in the chemokine 
concentration. 

Internalization can induce quasi-periodical alterations of the chemokine concentrations. Assume a 
population of B cells with full CXCR5 expression and high concentration of CXCL13. The receptor 
will be internalized on all cells using up the CXCL13 to some extent. The lowered CXCL13 concen- 
tration and CXCR5 levels promote B cells to enter a random walk migration modus and spread. The 
low uptake of CXCL13 by the B cells with internalized CXCR5 permits the CXCL13 concentration to 
rise again. The combination of CXCL13 uptake and switch of cell migration between chemotaxis and 
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random migration leads to a local synchronization of the receptor dynamics of these cells. The spread- 
ing of the B cells destabilizes the shape of the PLF and becomes significant when local populations of 
B cells synchronously spread out. 

The instability of the B cell aggregate leads to high B cell densities outside the FDC network. If 
the FDC network dynamics is switched on again this culminates in the generation of new FDC at the 
border of the PLF. The area covered by the B cells becomes more extended than expected from a 
densely packed set of B cells in the FDC network (see also [23]). In other words, the volume of the 
FDC network becomes bigger than the volume required by the number of B cells. Thus, the whole 
PLF becomes unstable in shape and follows the shifted 'center of mass' of the FDC network as a result 
of the CXCL13 concentration peaks. 

The positive feedback loop (CXCL13 stimulating higher LTai/32 levels in B cells, inducing more 
CXCL13 secreting FDC [14]) enhances this effect because B cells at low density can already induce 
FDC when their LTai/32 levels are high due to the CXCL13 stimulus. This is most pronounced when 
the peak LTai/32 levels of B cells are delayed compared to the peak CXCL13 stimulus because then 
the cells have a high probability to enter a FDC free area with enhanced surface LTai/^ levels. 

Note that in a previous study no instability was induced by internalization dynamics [23]. This 
is related to SIP chemotaxis towards the LE and to a wrong position of the LE in this study: The 
follicle formed around the LE position. 

In conclusion, internalization dynamics of CXCR5 induce unstable PLF. As such instabilities are 
not observed in nature, the present model predicts that receptor internalization is not a dominant 
process in PLF formation and maintenance. Even though it can not be excluded that internalization 
is active and only attenuated by further interactions, it is unlikely that we loose essential features of 
PLF formation when neglecting internalization dynamics in the following. However, the robustness of 
the results with respect to internalization dynamics will be considered. 

3.2 Including lymphatic endothelium dynamics 

LE dynamics might be considered to stabilize follicle shape and size by locating the exit points for 
B cells. However, LE dynamics turn out to destabilizes the follicle system independently of the 
internalization dynamics (Fig. [H E-H and Fig- E])- The reason for this instability is as follows: B cells 
that approach the LE - either by random motion or SIP-directed migration (Sec. 13.4]) - can stimulate 
FRC to become FDC if they occur at sufficient density such that FDC are generated right next to the 
exit spots. As FDC and LE are anti-correlated, i.e. the vicinity of LE is free of FDC, the generation 
of FDC leads to the degradation of the LE. The exit spots for B cells are pushed further away. The 
B cells need to approach the LE and 'follow' the distant exit spots. Due to the increased distance of 
the exit spots from the follicle center the B cells take longer to reach the spots thus increasing their 
transit time. Consequently the follicle is enlarged and the density of B cells close to the exit spots 
again reaches concentrations to induce new FDC. Thus a self-perpetuating process drives the border 
between FDC and LE leading to a constantly growing PLF size. 

Even though the model for the genesis of lymph endothelium is rather simplistic, it nevertheless cov- 
ers a reasonable idea of the dynamics of LE. The destabilizing effect of this basic model points towards 
an additional process that prevents the PLF from constant growth by inhibition of the generation of 
FDC at the right follicle size. 
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Figure 2: The dynamics of the efferent lymphatic destabilizes the follicle shape. The limiting factor 
for the follicle growth is mainly the influx of B cells as shown by the number of cells at each 
time point. In the case of low density LE there exist metastable states where the B cell 
number is almost constant while the number of FDC is growing until it reaches the LE and 
pushes it away which then leads to increased B cell numbers. The number of FDC has been 
multiplied by a factor 500 for reasons of clear visibility. 

3.3 Negative regulation of FDC generation 

A negative regulation mechanism that inhibits the generation of FDC can be realized in two distinct 
ways. First, the FDC induction is prevented when some factor is missing, i.e. FDC precursor activation 
gets lost. That implies that cells located in the T zone but not in the follicle secrete a substance that 
keeps the FDC progenitors in a state in which they are able to become FDC. The production and 
sensitivity for the signal has to be such that only many cells are able to give rise to a sufficient signal 
concentration. The fact that the T zone can be absent while separated PLF are observable in certain 
knockout experiments [9] makes this scenario rather unlikely. 

The second option is that the generation of FDC leads to the secretion of a substance X that 
inhibits their own generation counter-balancing the effect of LTai/32- Note that the negative regulation 
shall not affect existing FDC and not induce their dedifferentiation to FRC. Also the signal X has 
to be diffusible. A local inhibition mechanism might prevent the generation of FDC but does it 
independently of the follicle size. A signal sensitive to the follicle size needs to be accumulating and 
dispersing the information in the vicinity of the PLF. 

The simulations show that negative regulations of FDC induction, indeed, leads to stable follicles. 
Both aforementioned scenarios have similar effects (Fig.[TJ I and J). The PLF forms close to the HEV as 
expected from previous work where a small area containing LE had been positioned in a relatively large 
distance to the HEV [23] . In this configuration the PLF were forming around the LE and the HEV. In 
the present simulation the formation around the LE is no longer possible, because the homogeneous 
distribution of LE does not provide an aggregation center for the B cells to reach sufficient high 
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density for induction of FDC via a large LTai/32-stimulus. The lack of a B cell aggregation center in 
the present work is rather independent on SIP acting chemotactically on B cells or not. The broad 
distribution of SIP-producing LE can only provide shallow gradients not leading to B cell aggregation. 
The chemotactic activity of B cells in response to CCL21 can however lead to elongated follicles when 
the a CCL21 concentration peak is nearby the PLF center (Fig. [TJ J). 

The presence of one of these negative regulation mechanism has impact on the morphology of the 
PLF now exhibiting two zones. One B cell rich zone with FDC and without vessels, and a shell of B 
cells with low or vanishing FDC density but with LE. Within this shell B cells are exiting from the 
PLF and SLT. 

An additional stability test can be performed using a preformed FDC network that is placed in a 
certain distance from the HEV. The simulation of this configurations shows that the preformed FDC 
network rapidly extents to its final size and most importantly remains at its location (Fig. [H K). Thus, 
in contrast to simulations in [23], the B cells passing through the LE and the FRC network are not 
able to induce FDC before approaching in the PLF. 

3.4 Why SIP chemotaxis is not likely to occur in primary lymphoid follicles 

An important result of this study is that the chemotaxis towards SIP is no longer required and 
even more that SIP chemotaxis might not be active in the PLF system. In contrast to the previous 
study [23] the SIP chemotaxis is not required for B cells to find an exit spot because the LE is 
distributed homogeneously around the follicle. This makes it easy for B cells to find the LE by 
random migration within a rather short time. Thus, the role of a chemical signal has been replaced 
by a specific functional morphology. 

The effect of SIP in the PLF system is mainly a disturbed border of the lymphoid follicle (Fig. [TJ L). 
This results from a rather low concentration of CXCL13 at the boundary of the follicle which allows 
SIP to influence the direction of motion of B cells at the follicle border. This is only slightly affected 
by the internalization dynamics of CXCR5 because the cells at the border are mainly sensitive for 
CXCL13 due to low CXCL13 concentration. The morphology caused by SIP suggests that B cells in 
the follicle are not chemotactically responding to SIP. 

4 Discussion 

An agent-based lattice-free simulation tool for the formation and maintenance of PLF was presented. 
The generation of FDC is induced by B cells providing sufficient LT«i/?2 signals to stromal progenitor 
cells. B cells are sensitive to chemokine CXCL13 via their receptor CXCR5. CXCL13 is secreted by 
FDC. The positive feedback loop that enhance B cell LTai/^2 levels by CXCL13 stimulation [14] is 
not required to understand PLF formation. Optionally, the desensitization of the CXCR5 receptor 
via ligand- induced receptor internalization is considered. A simple model for the microanatomy of the 
lymphatic vessels in the vicinity of lymphoid follicles is used assuming an anti-correlation between LE 
and FDC. 

In contrast to previous simulations [23], which did not consider CXCR5 internalization dynamics 
or LE anti-correlation with FDC, PLF are found to be unstable when at least one of these additional 
features is included. In the following this is explored in more detail, and a mechanism is proposed to 
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elevate the problem of follicle instability. This model prediction is discussed in the context of related 
experiments. 

Internalization of chemokine receptors does not drive PLF formation 
The instability generated by the internalization dynamics of CXCR5 on B cells in the simulation 
suggests that this process might not be relevant in the PLF system. The internalization of CXCR5 
has been introduced in the simulation to desensitize chemotactic attraction of B cells to the center of 
the FDC network and, thus, to make them approach the follicular border either by random migration 
or chemotactic migration in response to other chemokines like CCL21, CCL19, and/or SIP, and, 
ultimately, to allow B cells to leave the follicle. A couple of points argue against this process being 
relevant to B cells leaving the follicle. 

In order to compensate for shape fluctuations induced by the CXCR5 internalization dynamics and 
achieve a spherical PLF the production of CXCL13 could in principle be raised. The resulting CXCL13 
concentration has to be high enough to ensure that CXCL13 removal by CXCR5 internalization is 
not strong enough to invert the CXCL13 gradient. However, the high CXCL13 concentration would 
result in high CXCL13 concentration everywhere in the PLF leading to rapidly CXCR5 desensitized 
B cells even at the border of the PLF. Consequently CXCL13 would not be able to hold B cells for a 
long time and B cells would leave the follicle after 15-30 minutes [67]. This short PLF transit time is 
in contradiction to the experiments observing transit of at least 3-5 hours [33,43-45]. 

Another argument against the internalization dynamics is that it cannot be stable in regulating the 
size of follicle. Assume that by random fluctuation more B cells arrive at the follicle. Internalization 
reduces the CXCL13 levels such that cells are less desensitized and remain longer in the follicle. As a 
consequence the follicle is growing in B cell numbers. The generation of FDC (with and without the 
positive feedback loop [14]) is by far slower than cell migration such that the B cells number remains 
at a higher ratio compared to the CXCL13 sources and letting the follicle grow further. Similar 
arguments apply to the case when a fluctuation leads to lower B cell numbers. When less B cells 
are in the follicle, then CXCL13 is dominating and B cells leave the follicle fast by a more rapid and 
frequent CXCR5 internalization. The CXCL13 concentration in the shrinking follicle is rising further 
until the slower dynamics of the FDC reduces the CXCL13 amount. However, due to their faster 
dynamics the number of B cells decreases faster. Basically, the slow dynamics of the FDC compared 
to the migration and internalization dynamics of B cells is responsible for this instability. 

Note, that it cannot be ruled out that a another modification of the CXCR5 receptor levels, for 
example on the transcriptional level, leads to a desensitization of B cells for CXCL13 thus enabling 
them to leave the PLF. Also it cannot be excluded that mechanisms exist that compensate for the 
internalization induced instability of the PLF. This may be a desired situation considering that the 
internalization dynamics can induce cycling of naive B cells between the center of the PLF and its 
border. The cycling is due to the strong desensitization in the center of the PLF where high CXCL13 
concentrations exist. B cells can then approach the border by random migration where the lower 
CXCL13 concentrations permit a resensitization for CXCL13 chemotaxis. These B cells migrate back 
to the center of the PLF. 

Follicular B cells are not responsive to CCL19, CCL21, and SIP 
The morphological data suggest that B cells in the PLF are not chemotactical responsive for SIP. A 
superposition of the two chemoattractants SIP and CXCL13 generates a blurry PLF border which is 
not observed experimentally. This is in line with the recent observation that SIP is not a chemotactic 
factor in vivo for T cells in the LN [27]. The same blurry PLF border would be induced if follicular 
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B cells responded to the T zone chemokines CCL19 and CCL21. Therefore, the simulation suggests 
that chemotaxis to CCL19, CCL21 and SIP is suppressed in follicular B cells, which seems to be 
the case in vivo [12,77,87,88] despite the positive response in vitro [74,76]. This may be caused by 
internalization of CCR7 at the border of the PLF where the concentrations of CCL19 and CCL21 are 
still quite high [82, 89] which is supported by the missing CCR7 surface levels of freshly isolated B 
cells [73]. 

A CXCL13-LTa±P2 positive feedback loop initiates PLF formation and gradual FDC differentiation 
The function of the positive feedback loop [14] (LTai/32 levels are increased by CXCL13, and LTai/?2 
induce novel FDC which are the source of CXCL13) was analyzed in the simulations. There are 
two possible functions: First, a randomly formed B cell aggregate may induce a minimum level of 
CXCL13 that starts the feedback loop leading to increased LTai/32 and finally to increased CXCL13. 
The CXCL13 levels reached in this way are sufficiently strong to induce chemotaxis in B cells and to 
dominate other chemotactic responses (e.g. to CCL21). This initiates formation of follicles. However, 
in vitro data show that very high CXCL13 levels (10 nM) are required to efficiently induce LTai/?2 on 
B cells [89]. These concentrations are already optimal for chemotaxis [90,91], which sets in question the 
former functional interpretation of the feedback loop. If one, nevertheless, assumes that the positive 
feedback loop is required to induce initial CXCL13 expression, then the initiation of the first FDC 
requires either another helper cell type or high LTai/?2 levels on some B cells. The latter might be a 
result of stimulation with antigen [14]. Then, either the cellular helper source needs to be identified 
or the dependence on B cell activation has to be confirmed. 

An alternative interpretation is that the LTai/?2 signal resulting solely from the aggregation of B 
cells may be sufficient to induce CXCL13-producing FDC but not other FDC markers like CD21 or 
CD35 which require the high LTai/?2 levels achieved by the feedback loop [14]. The gradual maturation 
of FDC observed in other studies [20,51] supports this notion. 

Moreover the notion of gradual maturation favors the strict interpretation of FDC as excited state 
of stromal cells and speaks against the differentiation of an unknown progenitor. This may also 
explain why in certain experiments FDC are rapidly lost when the LTai/^ stimulus is blocked with 
antibodies [13,62,92,93]. Thus the loss of the LTai/32 signal may not lead to the loss of FDC by 
apoptosis but a 'decay' of exited stromal cell back to their native state. 

Lymphatic vessels 

The microanatomical location of the LE has been modeled by an exclusion principle: Whenever a 
FDC is newly appearing the lymphatic vessels in the vicinity are removed. Presumably, remodeling 
is organized by the exchange of chemicals inducing behavioral changes of the LE cell. This is not 
considered explicitly but summarized in a time delay of LE remodeling of several hours. Analogously 
the dedifferentiation of a FDC into a FRC is followed by the creation of LE. 

The regulation of the LE is not precisely understood. In the following several hypothetic pathways 
that may regulate the LE dynamics and their corresponding experimental evidence are discussed. The 
LE dynamics could be related to the ECM. The reticular fibers are sparse in the follicle [31,94-100] 
and their degradation in the PLF could cause the disintegration of lymphatic vessels. A candidate 
signal to be involved is the expression of matrix metalloproteinases by FDC [101]. The generation of 
ECM components by FRC provides the ground for LE cells to extent the network in this area. The 
required vascular endothelial growth factors VEGF-C and/or VEGF-D can be provided from cells like 
macrophages (which can transdifferentiate into LE cells themselves) [102] or dendritic cells [103]. In 
this view, the gradual destruction of the reticular network in developing follicles in sheep [18] would 
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be accompanied by a gradual loss of LE in this area, which, in turn, promotes PLF growth. 

There exists indirect evidence how LE remodeling might be regulated in the PLF. Angiogenesis 
and lymphangiogenesis during airway infections is dependent on B cells [104]. A candidate for the 
mediator signal is vascular endothelial growth factor (VEGF)-D [104]. Although the cellular source 
is not clear VEGF-D could directly couple to B cells as suggested by the stimulation of VEGF by 
c-Myc + B cells in lymphoma [105]. A VEGF source candidate in SLT are DC in the T zone [103]. The 
co-localization of DC in the T zone correlates with the presence of LE. Thus it may turn out that DC 
are important for the structure of SLT. The lack of LE in the PLF may then be coupled to the absence 
of DC and not to a LE-destructing process initiated by the FDC. This hypothesis is already covered 
by the present simulations outcome of the simulation because in the real tissue as soon as FDC form, 
the DC are replaced by B cells. Thus VEGF will no longer be present and consequently lymphatic 
vessels are not maintained. An experiment using a system with VEGF-/- DC may elucidate the role 
of DC in SLT lymphangiogenesis. 

Assuming that DC are responsible for the location of LE, the ring of B cells found in certain LTai/?2 
and/or TNF-a-deficiencies [9-11] should be devoid of LE. However, the absence of LE in the B cell 
ring does not exclude the possibility that high B cells densities lead to the reduction of LE. Thus, the 
anti-correlation of FDC and LE assumed in the simulation can be re-interpreted to be indirect in the 
way that FDC support the aggregation of B cells leading to reduced LE in the PLF. 

TNF-a might be not only involved in the generation of FDC [9-11] but directly in the LE remodeling 
process. Both B and T cells can provide TNF-a [9,11]. TNF-a can upregulate VEGF-C ( [106] and 
references in [103]). The anti-correlation between FDC and vessels may be related to the TNF-a levels 
in the PLF that are raised by TNF-a-producing FDC [107]. It is known that TNF-a can also have 
anti-angiogenic effects at high concentrations [108] . TNF-a in the PLF may be present at sufficiently 
high levels to induce anti-angiogenic effects . 

The presented model assumes that the lymphatic network formation in SLT is coupled to the 
dynamics of PLF. In contrast one could assume that the LE is formed before PLF formation and 
determines the position and shape of the developing PLF. A direct test for the existence of the LE 
dynamics is to provide anti-angiogenic factors to LTai/3 2 or similar knockouts prior to reconstitution 
by wildtype lymphocytes or bone marrow. The PLF that form should be much smaller than usual and 
contain LE. Alternatively, the inverse experiment could be done by blocking LTai/?2 in wildtype mice 
or reconstitute them with bone marrow from LTai/3 2 • The disruption of the PLF should be followed 
by the presence of vessels in these areas. When using anti-angiogenic factors during the experiment, 
the gaps in the vessel network due to the PLF should be preserved when the PLF disappear. 

A novel mechanism of negative regulation of FDC generation 
The instabilities created by the LE dynamics require a mechanism that negatively regulates the induc- 
tion of FDC. Even though negative regulation might be realized in a number of different mechanisms, 
from the simulation viewpoint, only the source and the propagation of the signal initiating negative 
regulation are relevant. Either a signal has to be reduced or to be produced when FDC are generated. 
In both cases the signal overcomes a threshold at the border of the PLF. The signal is required to be 
diffusive in order to act on the length scale of a PLF and can be shown to require a rather fast decay 
in order to establish a time- independent (stable) follicle size (Appendix |A"|) . 

Within the framework of the simulation an inhibiting signal is suggested that is produced when 
FDC are generated. Note that this inhibiting signal may be secreted either by B cells located in the 
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PLF or by FDC. If the B cells are the sources of the hypothetic signal, they would constitutively 
produce the substance such that small B cell aggregates induce FDC while larger ones inhibit FDC 
induction. Alternatively, B cells sense the presence of FDC in the PLF (e.g via CXCL13) and produce 
the inhibitor upon stimulation in concurrence to LTctifo- 

The constitutive production of the inhibiting signal by B cells implies that the PLF would always 
grow from low B cell numbers to its final size and large aggregates that formed due to some other 
process would prohibit follicle formation. Considering the PLF formation in mice lacking T cells [9] 
this seems an unlikely situation because the B cell dense tissue clearly shows PLF formation. Also 
in normal MALT and LN organogenesis an intermediate ring of B cells around the developing T 
zone appears [6-10,109-112]. The ring formation suggests that high B cell numbers do not suppress 
FDC generation. Thus most likely either only the B cells in the follicle get stimulated to produce a 
negatively regulating factor or the FDC themselves are the source. 

Primary lymphoid follicles exhibit two zones 
The regulation of the follicle size under the influence of the dynamics of efferent lymphatic vessels 
implies that the B cell follicle consists of two zones. One zone contains B cells and the FDC network 
the other one B cells and efferent lymphatic vessels. From the corrosion casts this organization cannot 
be determined [35-37] but studies with FDC marker seem to support this view [10,14,20,92,112-115]. 
However, in the framework of the present simulation the FDC network is defined by the CXCL13 
distribution. Thus, the state of seemingly absent FDC in the outer zone of the PLF is not resolved in 
the model. The results may also indicate that the state of 'FDC-ness' decreases radially towards the 
border of the follicle. 

A theory of PLF formation and maintenance 
The simulation of PLF formation and maintenance gives rise to the following vantage point: FDC 
are a differentiated state of stromal precursor cells. Differentiation is promoted by LTai/^2, which is 
provided by B cells. The initiation of PLF formation is related to small aggregates of B cells and does 
not rely on the a positive feedback loop [14] between CXCL13 and LTai/32- B cells in a follicle have to 
desensitize their chemotactic responsiveness to T zone derived chemokines (CCL19, CCL21) in order 
to establish a realistic flow of B cells through the follicle and to avoid follicle instabilities. The model 
predicts that a factor negatively regulating FDC generation determines the size of a PLF. This factor 
stabilizes PLF shape and size. It is most likely an inhibitory diffusing signal secreted from either FDC 
or B cells residing in the PLF. 

A drawback of the presented model is that other desensitization mechanisms for the chemotactic 
response of B cells are not considered due to a lack of corresponding data. Other regulatory pathways 
like the transcription of the CXCR5 receptor may modify the PLF formation theory drawn in this 
study. The lymphangiogenesis model used here is very simple and only phenomenological. A detailed 
modeling of the lymphangiogenesis during PLF formation is left for future research, especially the use 
of a real vessel structure instead of spherical representatives. The dynamics of the adjacent T zone 
may influence PLF formation, although preliminary results seem not to change the conclusion drawn 
from the presented results. In addition to a more detailed lymphangiogenesis the dynamics of the 
HEV during PLF formation might be considered. 
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5 Method 



The PLF system is simulated using an agent-based off-lattice model [25,116-121]. To achieve a lattice- 
free description a regular triangulation has been chosen [23-26]. The regular triangulation, which is 
a generalization of the Delaunay triangulation, is used to provide the neighborhood topology for the 
cells that allows for a continuous representation of cell positions and sizes. The simulations of cells 
is represented in a 3-level model. The first level is the internal state of the cells representing the 
dynamics of the phenotype of the cell. The second level models the contact interaction between cells 
including mechanical interactions with the environment and exchange of signals by membrane bound 
molecules like LTai/?2. The third level incorporates long range interactions via diffusive substances 
which is used to describe the chemokine distribution in the PLF. The corresponding equations are 
briefly summarized in this section. More details of the model can be found in [23]. 

5.1 Internal cell dynamics 

The phenotype of a cell is described by a set of internal cell variables. These include internal times 
to indicate when which type of event may happen. An example is the persistence time of cells during 
chemotactic motion [75,82,85,86,122,123]. When the persistence time has past, the cell can reorient 
its migration direction to the local chemokine field [124,125]. Another variable is the differentiation 
of a FRC to a FDC where two times are important: The time the LTai/?2 stimulus is provided and 
the time the differentiation is completed (see below). 

A set of variables describes the mechanics of the cell, e.g. velocity, cell polarization, and cell volume. 
These variable couple directly to the next level of description, the contact interaction of cells, and 
have no influence on other internal states directly but can influence the variables of neighbor cells. 
Additionally there are parameters like the elasticity constant mentioned below that can be cell-type 
specific. Further details of the model and an extensive set of parameters have been published [23] and 
can be found in the supplementary material. 

Upon contact, two cells can exchange signals via their contact surface. In a simple approach with un- 
polarized ligand/receptor distributions the signal strength is proportional to the contact area provided 
by the Voronoi tessellation [23,25, 126]. 

5.1.1 Internal dynamics of FRC and FDC 

Cell differentiation is described as phenotype change of a cell. The realization of cell differentiation is 
derived from the change of cell behavior and internal states upon internal or external stimuli. Within 
the model a FRC differentiates to a FDC when the signal threshold for LTai/32 has been exceeded 
for a given time Tfrc^fdc- The signal is determined by summing up all LTai/32 contributions 
from neighbor cells, i.e. surface density of LTai/?2 times contact area. The differentiation is then 
instantly performed changing the internal cell states of a FRC into that of a FDC. In a similar 
manner FDC differentiate back to FDC after the LTai/?2 signal is below the threshold for a given 
time Tfdc^frc- T° reduce the amount of parameters and considering the lack of experimental data 
Tfrc<-fdc = ?fdc^frc is assumed. Also the thresholds for LTai/32 are identical constants for FRC 
to FDC and FDC to FRC differentiation. 

The positive feedback loop [14] involving B cells is realized by a linear interpolation of surface 
LTai/32 level between a constitutive low level and CXCL13-induced high LTai/32 level. The induction 
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by CXCL13 is chosen to be proportional to the fraction of bound CXCR5 receptor implying an 
immediate response. As an additional variant the LTai/32 level can been chosen proportional to the 
amount of internalized CXCR5 to model a delayed response. Additional delays which may result by 
the transcription of LTai/?2 have not been considered. This seems to be realistic considering that B 
cells reside in the PLF for a few hours [33, 43-45] which permits high LTai/?2 levels only when the 
cells are about to leave the follicle, if at all. 



5.2 Equations of motion 

The contact interaction of cells is primarily the mechanical interaction between cells. It is described by 
Newtonian equations of motion in the overdamped approximation. In this approximation acceleration 
of cells and consequently conservation of moment can be ignored. The equations incorporate elastic cell 
responses upon cell contact, forces generated by actively migrating cells, and friction forces between 
cells as well as cells and extracellular matrix. 
The Newtonian equations of motion read 
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The forces at position of the cell i depend on the internal state <f> of the cells for active forces F act 
and on the cell's position as well as on the position of all neighbor cells j £ Mi for passive forces 
ppass These forces are counter-balanced by the velocity-dependent drag forces F drag resulting in a 
ODE system of first order for the cell positions. 

The passive forces are the elastic described by the JKR-model [127]. It depends on the virtual cell 
overlap hij = ri + rj — dij where dij is the center distance of the cells and ri and rj are the cell radii. 
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with cell elasticity constants Ei and Ej, Poisson numbers i/j and Vj. 

The force acting on cell i by exerting active forces on a neighbor cell j reads 



Fff(&) = (HjPi sign([x* i - x*] • a) 
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with the cell orientation Ojj , cell surface contact point x*^ , constriction ring center x* , interaction area 
dij and the pressure p* exerted by cell i. Additionally a constant active force — F act is directly exerted 
by the cell on the ECM by the cell driving the cell in the direction of the cell axis Oj . 
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The drag force is linear in the velocity v given by 

F? ^ = -VmedTi [ 1 " ) Vi 

o,- (4) 

jeAfi 4 

with specific coefficients r/j, and /^medium- ^ ot is the total surface of a cell and A = ^ ■ djj is the 
surface in contact with other cells. 



5.3 Reaction-diffusion system of chemokines 

The chemotaxis of cells is described by coupling the direction of the active force of a cell to the local 
chemokine gradient. According to the observation that leukocytes tend to have a persistence time 
between subsequent orientation changes [75, 82, 85, 86] the gradient is sensed by the simulated cells 
periodically. The concentration of the chemokines CXCL13, and CCL21 are computed solving the 
time dependent diffusion equation. CCL19 is not explicitly calculated as it has a similar distribution 
like CCL21 [12], acts on the same receptor [74], and has a far lower concentration that CCL21 [89]. 
Note, that in inhomogeneous media like complex tissue a reaction-diffusion approach can become 
invalid and refinements might be necessary [128]. 

The receptor dynamics of the free (R) and internalized receptor (R*) are described by the following 
equations: 

R = -k on R c + k oS (R tot -R-R*) + k v R* 

R* = h(R tot - R - R*) - k T R* (5) 
c = -k on Rc + k s (i?tot - R - R*) - kc + Q. 

Rtat represents the constant number of receptors in any state of a cell. k on and k Q s are the rate 
constants for the ligand receptor binding. The variable c describes the concentration of the ligand 
which decays with the constant k and is produced by the source Q. k\ and k r are the rate constants 
for internalization of the receptor-ligand complexes and the recycling of receptors, respectively. 

In the model the LE is assumed to attract lymphocytes via the chemoattractant SIP. For simplicity 
and based on the notion that the dynamics of the LE is coupled to a slow process - the FDC dynamics 
- and by itself is a slow process it is calculated using the Poisson equation. Thus the time dependence 
of the SIP concentration is assumed to play no role and is omitted. Similarly the internalization 
dynamics of the SIP receptor is kept simple without explicit internalization dynamics as presented in 
(Eq. [5]). The reexpression of SIPi is made independent of the SIP levels and increases linearly until 
the levels permits the cell to enter LE after about 3 hr [33,42-45]. 



5.4 Initial conditions 

For most of the simulations a stromal background consisting of CCL21 producing FRC with a centrally 
localized preformed FDC network has been used. The area containing the FRC contains a constant 
density of exit spots which resemble the LE where cells can leave the tissue. A fixed number of HEV is 
located acentrically outside the preformed FDC network. The HEV are represented as spheres where 
cells enter the tissue. The majority of the simulations is performed with an influx of B cells only. 
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Table 1: The parameters used in the simulation of the PLF. The references given support the used 
value. If no reference or comment, the parameter is a systemic model parameter chosen to 
get sufficient accuracy in the simulation. 



19 



A Estimating the life time of a signal for the negative regulation of FDC 
generation 

Solving the stationary diffusion equation with a decay term and a homogeneous spherical source / of 
size L (spherical coordinates with radius p) 

= D {d^ + 2p d~p)~ KC ~ f@{L ~ p) (6) 
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with I = \JD I k the typical length scale of the system. The negative regulation should set on around 
the coordinate L where the source of the negative signal vanishes. The highest values at this border 
will be reached when L » I. However, the PLF growth will stop as soon as this threshold is reached 
such that the typical radius of the PLF will be in the order of a few I. For typical parameters of 
the diffusion coefficient of 10 . . . 100 fim 2 s _1 the value for k should be in the range of 10~ 3 . . . 10~ 2 s, 
i.e. the typical lifetime of the molecule mediating the negative regulation is in the order of minutes to 
tens of minutes. The decay term may be just the uptake of the signal by the cells it is acting on or 
its enzymatic inactivation. 
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